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CHARACTERISTICS OF LARGE SPACE SYSTEMS 


Proposed large space structures have many characteristics that make them diffi- 
cult to analyze and control. They are highly flexible - with components mathemati- 
cally modeled by partial differential equations or very large systems of ordinary 
differential equations. They have many resonant frequencies, typically low and 
closely spaced. Natural damping may be low and/or improperly modeled. Coupled with 
stringent operational requirements of orientation, shape control, and vibration 
suppression, and the inability to perform adequate ground testing, these character- 
istics present an unconventional identification and control design problem to the 
systems theorist. 

This presentation describes some of the research underway within Langley's 
Spacecraft Control Branch, Guidance and Control Division aimed at developing theory 
and algorithms to treat large space structures systems identification and control 
problems. The research areas to be considered are Computer-Aided Design Algorithms, 
and Systems Identification and Control of Distributed Systems. 


• Highly flexible 

• Many resonant frequencies 

• Low natural damping 

• Stringent operational requirements 

• Orientation 

• Shape 

• Vibration suppression 

• Limited ground testing 
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ORACLS UPDATE 


The established, tested computer-aided design system entitled "Optimal Regula- 
tor Algorithms for the Control of Linear Systems" (ORACLS) is being updated and 
modified so as to more easily accommodate the numerically difficult characteristics 
of large space systems lumped models. Modifications include greater use of UNPACK 
software (ref. 1) and inclusion of more robust Riccati equation algorithms (ref. 2). 
The ORACLS package (ref. 3) is also being expanded to allow multivariable frequency 
domain analysis and modern approaches to order reduction. 


OPTIMAL REGULATOR ALGORITHMS 
FOR THE CONTROL OF LINEAR SYSTEMS 

A modern control theory design package 
for time-invariant linear systems 

• Modular construction 

• Efficient numerical methods 

• Unified Continuous and discrete systems 

Constant and time - varying gains 
Deterministic and stochastic 

• Quadratic synthesis methods 

• Cosmic, NASA TP 1106, Marcel Dekker 
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FREQUENCY RESPONSE ANALYSIS 


A general purpose multivariable frequency domain analysis package (FREQ) has 
been constructed to be compatible with ORACLS. Given a multivariable unity gain 
feedback loop around a cascaded design system and dynamic compensator, the FREQ 
package has options to compute (1) a variety of transfer matrices, with loop broken 
at input or output, (2) singular values and vectors for the matrices, and (3) multi- 
variable Bode-like plots with maximum/minimum singular values plotted against 
frequency. Typical use of this package is in the analysis/design of compensators 
for spillover control. 
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ADVANCED CTLTIVARIABLE ALGORITHMS 


Further needs in large space structures control design have motivated other 
algorithm development* Algorithms for computing multivariable invariant zeros 
(ref. 4) and for treating controllability/observability are included since finite 
element models can have transmission zeros near closely spaced open-loop poles if 
sensors and actuators are not properly selected. High-order design models and/or 
high-order compensators typically occur and require order reduction techniques 
(ref. 5). It is anticipated that these and other algorithms will be collected into 
a new computer-aided design package (ORACLS II) motivated by the needs of large 
space structures controller design. 


• Multivariable zeros and relative controllability/observability 

• Applied to sensor/actuator selection 

• Multivariable frequency response package (FREQ) 

• Applied to spillover control 

• Algorithms based on Hankel-Norm Theory 

• Reduced order model/ compensator 

• Combine algorithms into new CAD package 

• ORACLS n 
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SYSTEMS IDENTIFICATION AND CONTROL OF DISTRIBUTED SYSTEMS 


Although there are currently many theoretical and numerical difficulties asso- 
ciated with control laws designed with partial differential equation models, it is 
felt that sensor/actuator technology and on-board computer capability will eventu- 
ally be improved to the point where distributed parameter methodology can be applied 
to large space structures. Anticipating these developments, the Spacecraft Control 
Branch is supporting research into systems identification and control theory and 
algorithms based on partial differential equation models. The approach currently 
being considered is to apply spline-based Galerkin projection approximation methods 
(ref. 6) and multivariable identif ication/control theory to models generic to large 
space structures. 
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HOOP— COLUMN APPLICATION 


A particular example of this research is given by the derivation of a parameter 
and state estimation procedure for distributed systems which was demonstrated with a 
model generic to the hoop-column antenna. The reflector surface of the antenna is 
assumed to be approximated by the static two-dimensional stretched membrane equation 
with appropriate boundary conditions and variable stiffness. The problem considered 
was to estimate the system state (u) and stiffness (E) from given applied force and 
displacement measurements • 


Static two-dimensional stretched membrane 


7 TF [rE<r.e>|^<r.ei] f ^ ,r. 8)1 = Nr. 6 ) 


Over Q = [e, R] x [o, 2 7r] 
with boundary conditions 

u(£,0)=Uq. u(r,o) = u ( r . 2 7r) 
u (R, 0) = 0 

• Given f(r,0 )and measured displacements 

u m ( r i’ e j)’ (i = 1 . L r ; i = 1 • • • * l e ] ’ estimate 

E (r, 0)and u (r, 0) within Q 
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SYSTEM IDENTIFICATION APPROACH 


The general approach for identification can be described as follows. After 
formulating a distributed model for the dynamic system, some mathematical approxi- 
mation technique, such as finite elements or splines, is used to project the 
identification problem onto a finite-dimensional subspace. The finite-dimensional 
identification problem is then solved within the subspace. The approximation is 
successfully refined and the identification problem repeatedly solved to produce a 
sequence of estimates to be analysed for convergence. This approach, when 
specialized to the hoop-column application, is outlined below. Details of the study 
may be found in reference 7. 


General 

• Distributed parameter formulation 

• Project onto finite dimensional 
subspace 


• Solve identification problem 
within subspace 

• Successively increase subspace 
dimension and solve identification 
problem to generate sequence of 
estimates 


Hoop-column application 

• Stretched membrane equation 

• Galerkin projection with linear 
spline basis functions 
Parameterize E(r,0) via cubic 
splines 

• Output error identification 
technique 

• Numerical algorithm 
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OPTIMAL LINEAR REGULATOR ALGORITHMS 


The same general approach is being investigated for controller design. When 
applied to dynamic systems, the spline-based Galerkin approximation procedure pro- 
duces high-order state equations. If linear quadratic regulator theory is applied 
for control law design, the related difficulty of solving high-order Riccati equa- 
tions arises. In control applications, wherein only the regulator gain is required, 
direct application of the Chandrasekhar technique (ref, 8) is extremely difficult due 
to numerical stiffness problems. Recent results from a NASA-sponsored research 
grant (NAG— 1—517) with Brown University have produced a new hybrid Chandrasekhar-type 
algorithm for approximating the steady-state gain matrix independently of the 
Riccati equations. 


• Linear quadratic regulator problem 

x = Ax+ Bu 

OO 

J = f (x l Cx + u T u) dt 
0 

A T P + PA - PBB T P + C T C =0 

u = -Kx = -b'px 

• Chandrasekhar algorithm 

-4 P(t) = A T P(t) + P(t)A 

- P(t)BB T P(t) + C T C, P(t f ) = 0 
-jp = -B T L T (t)L(t), K(t f ) = 0 

-L(t) [A - BK(t)] , L(t f ) = C 
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HYBRID ALGORITHM 


In order to begin the hybrid algorithm, the standard Chandrasekhar algorithm is 
employed to obtain an initial stabilizing gain for application of the Newton-Klein- 
man algorithm for solving the algebraic Riccati equation. After several steps 
through the Newton-Kleinman sequence, enough data are obtained to rewrite the recur- 
sion equation in an alternate form giving an update formula for the gain matrix. 
Banks and Ito at Brown University have discovered a way to compute the gain update 
(K^ + j) without computing Z^. 


Chandrasekhar equations give starting value 
for Newton-Kleinman 

Newton-Kleinman 

(A-BK-) T Pj + Pj (A-BKj) + kI Kj + C T C 


= 0 


K i + l =BP i 


Alternate form - 


(A-BKj) T Zj + Zj (A-BKj) + (Kj -K, _j_> T (K j -Kj^) 


= 0 


Kj+1 = Kj + B Z j 

Can compute update without computing Z, 
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HYBRID ALGORITHM (CONT'D) 


Apply the (Smith) bilinear transformation to the equation. A sequential 
solution can then be written for the tranformed equation from which the update 
algorithm can be derived. 


Bilinear transformation (Smith) - 


x = z i 

A = A - BK. 

D = K i ' K i-1 


U= ( I -rA ) _1 ( 1+rA) 

Y = 2r ( I -rA ) T D T D ( I -rA ) -1 


X = U XU + Y 


Sequential solution - 

x w = u V + v 

x k + i - x k = ° T( x k- x k-i> u 
Algorithm - 

X k+1 = X k +2rM k+l M k+l 
M k + 1 = M k U 

BTx k-i =B V 2rB HIi M k + i 
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NUMERICAL RESULTS 


For a distributed parameter example, consider the parabolic system below with 
boundary control. The control input u(t) is to be chosen to minimize J(u) subject 
to the dynamic equation and boundary condition constraints. After discretizing with 
linear splines and Galerkin projection, an Nth order linear quadratic regulator 
problem is obtained. The index N increases with refinement of the linear spline 
approximation. For N=10, CPU times are comparable between the Potter (ref. 8) and 
hybrid methods. However, as N increases, the hybrid algorithm excels. 
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N = 10 

0.162 

0.197 

N = 40 

6.3 

1.29 
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STABILITY AUGMENTATION BY BOUNDARY-FEEDBACK CONTROL 


Another approach to distributed-parameter formulation of the control problem is 
being considered by Balekriebnen (ref. 9) for application to the Spacecwf? SonSoJ 
laboratory Experiment (SCOLE) configuration (ref. 10). An abstract wave formulation 
is developed as a nonlinear wave equation in Hilbert space. The system is shown to 
be controllable, and a feedback control law is developed assuming point actuators 
and sensors at the boundaries. The control law is shown to be strongly stable and 
robust to parameter uncertainties. 


• M x(t) + A x(t) + K (x(t) ) + B u ( t ) = 0 

M: mass matrix 

A ; D-^L 2 [o,l] 3 xR 14 , differential operator 

K: nonlinear function 

u(t): applied moments, proof mass forces 

X= [V J * V- 1 ' V-’’ % (on u a * 0+, > u fl d-), 

u;(0+). u'g <0 + ), (0 + ), u;(L-), u' e (L-), u^fL-l, 

u <t* ( s 2 )’ u a ( s 2 )■ u <i> ( s 3 ) * u a( s 3)] T 
u ip < s )» u 0 (s): displacements; u (s); rotation 
s 2 > s 3 : P ro °i mass actuator locations 

• State space form: y(t) = Ay(t) + Kly(t)) + B ultl 

• Control law: u (t) = -P B*y(t) ; P>0 
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damping models 


The mathematical modeling of damping mechanisms is a little understood, most 
important aspect of systems identification of large space structures. Tradition- 
ally, damping is superimposed in linear viscous form on modal dynamic equations 
obtained from finite element software. Unfortunately, none of the expected forms of 
energy dissipation is linear viscous (ref. 11). Research employing partial 
differential equation dynamic models indicates that damping operators should have 

nonlocal structure. 


• Little understood 

• Recognized as most important aspect 

• Linear-viscous 

• • • ^ 

Mu + Du + Ku = f 

• None of expected damping is linear-viscous 

• PDE damping operators can require nonlocal structure 
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RESEARCH IN MODELING OF DAMPING 


Using the torsional beam equation to model the Spacecraft Control Laboratory 
Experiment (SCOLE) (ref. 10), it has been found that, for proportional damping to 

consideration of the total beam length in the damping operator formulation is 
required. Future research under the Brown University grant (NAG-1-517) will con- 
sider other forms of damping operators for an Euler-Bernoulli beam model. The 
investigation will include damping with nonlocal interaction terms and fading memory 
formulations where, for example, stress is a functional of the history of deforma- 
tion and prior times* 


Scale model ( Balakrishnan, UCLA) 
for infinite beam with torsion equation 
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Investigating (Banks, Brown Univ.) 
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In Euler-Bernoulli beam model 
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